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We review tlie recently found large-scale anomalies in the maps of temperature anisotropies in the 
cosmic microwave background. These include alignments of the largest modes of CMB anisotropy 
with each other and with geometry and direction of motion of the Solar System, and the unusually 

low power at these largest scales. We discuss these findings in relation to expectation from stan- 
dard inflationary cosmology, their statistical significance, the tools to study them, and the various 
attempts to explain them. 



I. INTRODUCTION: WHY LARGE SCALES 
ARE INTERESTING? 

The Copernican principle states that the Earth does 
not occupy a special place in the universe and that ob- 
servations made from Earth can be taken to be broadly 
characteristic of what would be seen from any other point 
in the universe at the same epoch. The microwave sky is 
isotropic, apart from a Doppler dipole and a microwave 
foreground from the Milky Way. Together with the 
Copernican principle and some technical assumptions, 
an oft-inferred consequence is the so-called cosmological 
principle. It states that the distributions of matter and 
light in the Universe are homogeneous and isotropic at 
any epoch and thus also defines what we mean by cosmic 
time. 

This set of assumptions is a crucial, implicit ingredi- 
ent in obtaining most important resiilts in quantitative 
cosmology. For example, it allows us to treat cosmic mi- 
crowave background (CMB) temperature fluctuations in 
different directions on the sky as multiple probes of a sin- 
gle statistical ensemble, leading to the precision determi- 
nations of cosmological parameters that we have today. 

Although we have some observational evidence that ho- 
mogeneity and isotropy are reasonably good approxima- 
tions to reality, neither of these are actual logical conse- 
quences of the Copernican principle. For example the ge- 
ometry of space could be homogeneous but anisotropic — 
like the surface of a sharp mountain ridge, with a gentle 
path ahead but the ground dropping steeply away to the 
sides. Indeed, three-dimensional space admits not just 
the three well known homogeneous isotropic geometries 
(Euclidean, spherical and hyperbolic - , and H^), 
but five others which are homogeneous but anisotropic. 
The two simplest are 5^ x and x E^ . These spaces 
support the cosmological principle but have preferred di- 
rections. 

Similarly, although the Earth might not occupy a privi- 
leged place in the universe, it is not necessarily true that 
all points of observation are equivalent. For example, 
the topology of space may not be simply-connected 
we could live in a three-dimensional generalization of a 



torus so that if you travel far enough in certain direc- 
tions you come back to where you started. While such 
three-spaces generically admit locally homogeneous and 
isotropic geometries, certain directions or points might be 
singled out when non-local measurements are considered. 
For example the length of the shortest closed non-trivial 
geodesic through a point depends on the location of that 
point within the fundamental domain. Similarly, the in- 
homogcneity and anisotropy of eigcnmodcs of differen- 
tial operators on such spaces are likely to translate into 
statistically inhomogeneous and anisotropic large scale 
structure, in the manner of Chladni figures on vibrating 
plates. 

The existence of non-trivial cosmic topology and of 
anisotropic geometry are questions that can only be an- 
swered observationally. In this regard, it is worth not- 
ing that our record at predicting the gross properties 
of the universe on large scales from first principles has 
been rather poor. According to the standard concordance 
model of cosmology, over 95% of the energy content of 
the universe is extraordinary dark matter or dark en- 
ergy whose existence has been inferred from the failure 
of the Standard Model of particle physics plus General 
Relativity to describe the behavior of astrophysical sys- 
tems larger than a stellar cluster — while the very homo- 
geneity and isotropy (and inhomogeneity) of the universe 
owe to the infiuence of an infiaton field whose particle- 
physics-identity is completely mysterious even after three 
decades of theorizing. 

The stakes are set even higher with the recent dis- 
covery of dark energy that makes the universe undergo 
accelerated expansion. It is known that dark energy can 
affect the largest scales of the universe — for example, 
the clustering scale of dark energy may be about the 
horizon size today. Similarly, inflationary models can in- 
duce observable effects on the largest scales via either ex- 
plicit or spontaneous violations of statistical isotropy. It 
is reasonable to suggest that statistical isotropy and ho- 
mogeneity should be substantiated observationally, not 
just assumed. More generally, testing the cosmological 
principle should be one of the key goals of modern obser- 
vational cosmology. 
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Section [V] categorizes and covers the variety of possible 
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explanations for these anomalies. We conclude in Sec. VI 



FIG. 1. The comoving length within the context of the concor- 
dance model of an arc seen at an fixed angle and the comoving 
Hubble length as functions of redshift. Linear perturbation 
theory is expected to work well outside the shaded region, in 
which the large scale structure (LSS) forms. 



With the advent of high signal-to-noise maps of the 
cosmic microwave background anisotropics and with the 
conduct of nearly-full-sky deep galaxy surveys, statisti- 
cal isotropy has begun to be precisely tested. Extraordi- 
nary full-sky temperature maps produced by the Wilkin- 
son Microwave Anisotropy Probe (WMAP), in particu- 
lar, are revolutionizing our ability to probe the universe 
on its largest scales [T}]S| . In the near future, these will be 
joined by higher resolution temperature maps and high- 
resolution polarization maps and, eventually, by deep all- 
sky surveys, and perhaps by tomographic 21-cm line ob- 
servations that will extend our detailed knowledge of the 
universe's background geometry and fluctuations into the 
interior of the sphere of last scattering. 

In this brief review, we describe the large-scale anoma- 
lies in the CMB data, some of which were first reported 
on by the Cosmic Background Explorer (COBE) Differ- 
ential Microwave Radiometer (DMR) collaboration in the 
mid 1990's. In particular, we report on alignments of the 
largest modes of CMB anisotropy with each other, and 
with geometry and direction of motion of the Solar Sys- 
tem, as well as on unusually low angular correlations at 
the largest angular scales. We discuss these findings and, 
as this is not meant to be a comprehensive review and we 
emphasize results based on our own work in the area, we 
refer the reader to literature for all developments in the 
field. This review extends an earlier review on the sub- 
ject by (author?) 7 , and complements another review 
on statistical isotropy in this Special Issue [5j. 

The paper is organized as follows. In Sec. |ll] we de- 
scribe the statistical quantities that describe the CMB, 
and the expectations for their values in the currently fa- 
vored ACDM cosmological model. 



In Sec. Ill we describe 
the alignments at the largest scales, as well as multipole 
vectors, which is a tool to study them. In Sec.|IV[ we de- 



ll. 



EXPECTATIONS FROM COSMOLOGICAL 
INFLATION 



scribe findings of low power at largest scales in the CMB. 



A fixed angular scale on the sky probes the physics of 
the universe at a range of physical distances correspond- 
ing to the range of observable redshifts. This is illustrated 
in Fig. [TJ where the comoving lengths of arcs at fixed an- 
gle are shown as a function of redshift, together with the 
comoving Hubble scale. Angles of 1 degree and less probe 
events that were in causal contact at all epochs between 
the redshift of decoupling and today; this redshift range 
includes physical processes such as the secondary CMB 
anisotropics. The situation is different for angles > 60 
degrees, which subtend arcs that enter our Hubble patch 
only at z < 1. Therefore, the primordial CMB signal on 
such large angular scales could only be modified by the 
physics of local foregrounds and cosmology in the rela- 
tively recent past (z < 1). Because they correspond to 
such large physical scales, the largest observable angular 
scales provide the most direct probe of the primordial 
fluctuations — whether generated during the epoch of 
cosmological inflation, or preceding it. 



A. Statistical isotropy 

What do we expect for the large angular scales of the 
CMB? A crucial ingredient of cosmology's concordance 
model is cosmological inflation — a period of acceler- 
ating cosmic expansion in the early universe. If we as- 
sume that inflationary expansion persisted for sufficiently 
many e-folds, then we expect to live in a homogeneous 
and isotropic universe within a domain larger than our 
Hubble volume. This homogeneity and isotropy will not 
be exact, but should characterize both the background 
and the statistical distributions of matter and metric fluc- 
tuations around that background. These fluctuations are 
made visible as anisotropics of the CMB temperature and 
polarization, which are expected to inherit the underly- 
ing statistical isotropy. The temperature T seen in direc- 
tion e is predicted to be described by a Gaussian random 
fleld on the sky (i.e. the 2-sphere S^), which implies that 
we can expand it in terms of spherical harmonics Ygrnis) 
multiplied by independent Gaussian random coefficients 
aim of zero mean. 

Statistical isotropy implies that the expectation values 
of all ri-point correlation functions (of the temperature 
or polarization) are invariant under arbitrary rotations 
of the sky. As a consequence the expectation of the tem- 
perature coefficients is zero, (aim) — 0, for all £ > 
and m = —£, ~£ + 1, . . . , +£. The two-point correlation 
becomes a function of cos = ei ■ €2 only, and can be 



expanded in terms of Legendre polynomials 
(r(ei)T(e2)) = C{e) ^ ^Y.^2t+l)CiPpXcose). (1) 



angular 2-point correlation 



Statistical independence implies that expectations of a^,„ 
with different I and m vanish. In particular, the two- 
point correlation function is diagonal in I and m 



(2) 



Statistical isotropy adds that the constant of proportion- 
ality depends only on ^, not m 



Im'^i'm') — 6ui5mm'Cl. 



(3) 



The variance Ci is called the angular power of the mul- 
tipole i. The higher n-point correlation functions are 
constrained in similar ways but, as we will see below, 
are not expected to provide independent information if a 
simple inflationary scenario was realized by nature. 



B. Gaussianity 

If inflation is driven by a single dynamically relevant 
degree of freedom with appropriate properties (minimal 
coupling, Minkowski vacuum in UV limit, etc.), then we 
can reduce the quantization of matter and space-time 
fluctuations during inflation to the problem of quantizing 
free scalar fields. For free fields the only non-trivial ob- 
ject is the two-point correlation (the propagator), and all 
higher correlation functions either vanish or are just some 
trivial combination of the two-point function. This prop- 
erty is mapped onto the temperature field of the CMB. A 
classical random field with these properties is a Gaussian 
with mean Tq and variance C{0). Thus the brightness 
of the primordial CMB sky is completely characterized 
by To and C{9) (or Ci). Note that evolution of pertur- 
bations leads to deviations from Gaussianity that would 
mostly be evident at very small scales {i ^ 100). More- 
over, many inflationary models predict small deviations 
from Gaussianity; these are described in other contribu- 
tions to this volume [HI HO] ■ 



Scale-invariance 




FIG. 2. Mean and (cosmic) variance of the angular two-point 
correlation function as expected from cosmological inflation 
(arbitrary normalization). Only statistical isotropy, Gaus- 
sianity and scale-invariance are assumed. Tensors, spectral 
tilt, reionization and the integrated Sachs- Wolfe effect are ne- 
glected for the purpose of this plot. Comparison to the pre- 
diction from the best-fit ACDM model (Fig. [5| reveals that 
these corrections are subdominant. Note that cosmic variance 
errors at different values of 6 are very highly correlated. 



of the Hubble parameter during inflation, which is mea- 
sured by the so-called first slow-roll function ei = du 
where dn = is the Hubble distance. From the weak 
energy condition ei > 0, while ei <C 1 during slow roll 
inflation. 

At the level of the angular power spectrum, exact scale- 
invariance implies the Sachs- Wolfe "plateau" (i.e. con- 
stancy of l{l + l)Ci at low t) [TT] 



2ttA 



(4) 



Here, again in the slow-roll parametrization, A ^ 
{Hinfi/Mp)'^T^ / ei. This neglects secondary anisotropics 
like the late time, integrated Sachs- Wolfe effect (particu- 
larly important at very low £) and the contribution from 
gravitational waves. Furthermore, inflation predicts a 
small departure from scale invariance, which has recently 
been detected (e.g. [5]), and which also contributes to a 
tilt in the aforementioned plateau. 



Another generic feature of inflation is the almost scale- 
invariance of the power spectrum of fluctuations. This 
can be understood easily, as the Hubble scale is ap- 
proximately constant during inflation as the wavelengths 
of observable modes are redshifted beyond the horizon. 
Given that fluctuations of modes on horizon exit are re- 
lated to the Hubble parameter, 5(j) = H/2tt, these modes 
have similar amplitudes. However, scale-invariance is not 
exact. In canonical slow-roll inflation models, the devi- 
ation from exact scale-invariance is due to the evolution 



D. Cosmic variance 

As we can measure only one sky, it is important to flnd 
the best estimators of Ci and C{0). Let us for the mo- 
ment assume that we are able to measure the primordial 
CMB of the full sky, without any instrumental noise. We 
also restrict ourselves to i' > 2, as the variance of the 
monopole cannot be defined and the measured dipole is 
dominated by our motion through the universe, rather 
than by primordial physics. (Separation of the Doppler 
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dipole from the intrinsic dipole is possible in principle 
[12j[T3], but not with existing data.) Statistical isotropy 
suggests to estimate the angular power by 



L=2 



L=3 



L=4 



L=5 



C, = 



2£ 



1 +^ 

— y 
+1 ^ 



(5) 



which satisfies (Cf) = Ci and is thus unbiased. The vari- 
ance of this estimator can be calculated assuming Gaus- 
sianity: 



Var(Q) 



2£+l 



-C 



(6) 



It can be shown that, assuming statistical isotropy and 
Gaussianity, Cg is the best estimator in the sense that 
it has minimal variance and is unbiased. However, we 
emphasize that these qualities depend intrinsically on the 
correctness of the underlying assumptions. 

With these same assumptions, the variance of the two- 
point correlation function is easily shown to be 



Var[C(6l)] 



1 

8^ 



^(2^-Hl)CjP/(cos0), 



(7) 



where C{0) is calculated from C( following Eq. ([ij). 

Putting the results of this section together allows us 
to come up with a generic prediction of inflationary cos- 
mology for C{9) on the largest angular scales; see Fig. [2] 



III. ALIGNMENTS 

In brief, the upshot of the previous section is that the 
twin assumptions of statistical isotropy and Gaussianity 
are the starting point of any CMB analysis. The mea- 
surements of the CMB monopole, dipole and (Ar)rnis 
tell us that isotropy is observationally established at the 
per cent level without any cosmological assumption, and 
at a level 10~^ if we attribute the dominant contribution 
to the dipole to our peculiar motion. For the purpose of 
cosmological parameter estimation, the task is to test the 
statistical isotropy of the CMB brightness fluctuations. 
At the largest angular scales, this can only be done by 
means of full sky maps. 

Let us assume that the various methods that have been 
developed to get rid of the Galactic foreground in single 
frequency band maps of the microwave sky are reliable 
(though we argue below that this might not be the case) . 
Our review of alignments will be based on the internal 
linear combination (ILC) map produced by the WMAP 
team, which is based on a minimal variance combina- 
tion of the WMAP frequency bands. The weights for the 
flve frequency band maps are adjusted in 12 regions of 
the sky, one region lying outside the Milky Way and 11 
regions along the Galactic plane. 
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FIG. 3. Multipole vectors of our sky, based on WMAP five 
year full-sky ILC map and with galactic plane coinciding with 
the plane of the page. The temperature pattern at each mul- 
tipole ^ (2 < ^ < 8) can either be described by an angular 
temperature pattern (colored lobes in this figure) , or alterna- 
tively by precisely I multipole vectors (black "sticks"). While 
the multipole vectors contain all information about the di- 
rectionality of the CMB temperature pattern, they are not 
simply related to the hot and cold spots and, for example, do 
not correspond to the temperature minima/maxima [14]. No- 
tice that 1 = 2 and 3 temperature patterns are rather planar 
with the same plane, and that their vectors lie approximately 
in this plane. Adopted from Ref. |15| . 



A. Multipole vectors 

To study the orientation and alignment of CMB mul- 
tipoles, Copi, Huterer & Starkman 115^ introduced to 
cosmology the use of multipole vectors; an alternative 
representation of data on the sphere. The multipole vec- 
tors contain information about the "directions" associ- 
ated with each multipole. In this new basis the temper- 
ature fluctuation multipole, may be expanded as 



rn — — ^ 



(8) 



where ti(^'*) is the i**^ vector for the multipole, e is 
the usual radial unit vector, Te is the trace of the pre- 
ceding product of multipole vector terms, and A^^"^ is the 
"power" in the multipole. By construction, we immedi- 
ately see that Ti is a symmetric, traceless, rank £ tensor. 
Subtracting 7} ensures that Ti is traceless and the dot 
products explicitly show this is a rotationally invariant 
quantity (a scalar under rotations). This form makes 
the symmetry properties obvious. As an example the 
quadrupole is written as 



(^)(24).e)(f)(2,2).e) 



:^)(2,i) . ^(2,2) 



(9) 



These two forms for representing T^, harmonic and 
multipole-vector, both contain the same information. At 
the same time, they are fundamentally different from 
each other. Each unit vector -O^^'*) has two degrees of 
freedom while the scalar, ^^^^ has one; thus the multipole 
vector representation contains the full 2£ + 1 degrees of 
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freedom. Note that we call each £i(^'') a multipole vector 
but it is only defined up to a sign. We can always reflect 
the vector through the origin by absorbing a negative 
sign into the scalar A^^^; thus these vectors actually are 
headless. Regardless, we will continue to refer to them as 
multipole vectors and not use the overall sign of the vec- 
tor in our analysis. This issue is equivalent to choosing 
a phase convention, such as the Condon-Shortley phase 
for the spherical harmonics. For the work reviewed here 
the overall phase is not relevant and thus will not be 
specified. 

An efficient algorithm to compute the multipole vec- 
tors for low-£ has been presented in [15] and is publicly 
available [TB]; other algorithms have been proposed as 
well p7HT5] . Interestingly, after the publication of the 
CHS paper [15], Weeks [18] pointed out that multipole 
vectors were actually first used by Maxwell [50] more 
than 100 years ago in his study of multipole moments in 
electrodynamics. They remain in use in geometrology, 
nuclear physics, and other fields. 

The relation between multipole vectors and the usual 
harmonic basis is very much the same as that between 
Cartesian and spherical coordinates of standard geome- 
try: both are complete bases, but specific problems are 
much more easily addressed in one basis than the other. 
In particular, we and others have found that multipole 
vectors are particularly well suited for tests of planarity 
and alignment of the CMB anisotropy pattern. More- 
over, a number of interesting theoretical results have been 
found; for example, Dennis [21 analytically computed 
the two-point correlation function of multipole vectors 
for a Gaussian random, isotropic underlying field. Nu- 
merous quantities have been proposed for assigning di- 
rections to multipoles and statistics on these quantities 
have been studied. In Copi et al. |14j we have summa- 
rized these attempts and have shown their connections 
to the multipole vectors. 



B. Planarity and Alignments 

Tegmark et al. ^221 and de Oliveira-Costa et al. [23] 
first argued that the octopole is planar and that the 
quadrupole and octopole planes are aligned. In Schwarz 
et al. [24 , followed up by Copi et al. [H] [5S] , we inves- 
tigated the quadrupole-octopole shape and orientation 
using the multipole vectors. The quadrupole is fully de- 
scribed by two multipole vectors, which define a plane. 
This plane can be described by the "oriented area" vector 



(Note that the oriented area vector does not fully char- 
acterize the quadrupole, as pairs of quadrupole multipole 
vectors related by a rotation about the oriented area vec- 
tor lead to the same oriented area vector.) The octopole 
is defined by three multipole vectors which determine 
(but again are not fully determined by) three area vec- 




FIG. 4. Quadrupole and octopole {£ — 2 and 3) temperature 
anisotropy of the WMAP sky map in galactic coordinates, 
shown with the ecliptic plane and the cosmological dipole. In- 
cluded are the multipole vectors (solid diamonds); two for the 
quadrupole (red diamonds) and three for the octopole (green 
diamonds). We also show the four normals (solid squares) to 
the planes defined by vectors that describe the quadrupole 
and octopole temperature anisotropy; one normal is defined 
by the quadrupole (red square) and three by the octopole 
(green squares). Note that three out of four normals lie very 
close to the dipole direction. The probability of this alignment 
being accidental is about one part in a thousand. Moreover, 
the ecliptic plane traces out a locus of zero of the combined 
quadrupole and octopole over a broad swath of the sky — 
neatly separating a hot spot in the northern sky from a cold 
spot in the south. These apparent correlations with the solar 
system geometry are puzzling and currently unexplained. 

tors. Hence there are a total of four planes determined 
by the quadrupole and octopole. 

In Copi et al. [25] we found that (see Fig. [4| 

• the four area vectors of the quadrupole and oc- 
topole are mutually close (i.e. the quadrupole and 
octopole planes are aligned) at the 99.6% C.L.; 

• the quadrupole and octopole planes are orthogonal 
to the ecliptic at the 95.9% C.L.; this alignment 
was at 98.5% C.L. in our analysis of the WMAP 1 
year maps. The reduction of alignment was due to 
WMAP's adaption of a new radiometer gain model 
for the 3 year data analysis, that took seasonal vari- 
ations of the receiver box temperature into account 
— a systematic that is indeed correlated with the 
ecliptic plane. We regard that as clear evidence 
that multipole vectors are a sensitive probe of align- 
ments; 

• the normals to these four planes are aligned with 
the direction of the cosmological dipole (and with 
the equinoxes) at a level inconsistent with Gaussian 
random, statistically isotropic skies at 99.7% C.L.; 

• the ecliptic threads between a hot and a cold spot 
of the combined quadrupole and octopole map, fol- 
lowing a node line across about 1/3 of the sky and 
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separating the three strong extrema from the three 
weak extrema of the map; this is unHkely at about 
the 95% C.L. 

These numbers refer to the WMAP ILC map from three 
years of data; other maps give similar results. Moreover, 
correction for the kinematic quadrupole - slight modifi- 
cation of the quadrupole due to our motion through the 
CMB rest frame - must be made and increases signifi- 
cance of the alignments. See Table 3 of Ref. ^5, fo^' the 
illustration of both of these points. 

While not all of these alignments are statistically in- 
dependent, their combined statistical significance is cer- 
tainly greater than their individual significances. For 
example, given their mutual alignments, the conditional 
probability of the four normals lying so close to the eclip- 
tic, is less than 2%; the combined probability of the four 
normals being both so aligned with each other and so 
close to the ecliptic is less than 0.4% x 2% = 0.008%. 
These are therefore clearly surprising, highly statistically 
significant anomalies — unexpected in the standard in- 
fiationary theory and the accepted cosmological model. 

Particularly puzzling are the alignments with solar sys- 
tem features. CMB anisotropy should clearly not be cor- 
related with our local habitat. While the observed cor- 
relations seem to hint that there is contamination by a 
foreground or perhaps by the scanning strategy of the 
telescope, closer inspection reveals that there is no obvi- 
ous way to explain the observed correlations. Moreover, 
if their explanation is that they are a foreground, then 
that will likely exacerbate other anomalies that we will 
discuss in section ITV Bl below. 

Our studies (see [Mj) indicate that the observed align- 
ments are with the ecliptic plane, with the equinox 
or with the CMB dipole, and not with the Galactic 
plane: the alignments of the quadrupole and octopole 
planes with the equinox/ecliptic/dipole directions are 
much more significant than those for the Galactic plane. 
Moreover, it is remarkably curious that it is precisely 
the ecliptic alignment that has been found on somewhat 
smaller scales using the power spectrum analyses of sta- 
tistical isotropy 

Finally, it is important to make sure that the results are 
unbiased by unfairly chosen statistics. We have studied 
this issue extensively in |14j , and here we briefly describe 
the principal statistics used to quantify the probability 
of alignments quoted just above. 

To define statistics we first compute the three dot- 
products between the quadrupole area vector and the 
three octopole area vectors. 



A, = |u;(^^l'2).^(3;i,,)|^ 



(11) 



The absolute value is included since the multipole vectors 
are headless; thus each Ak lies in the interval [0, 1]. Two 
natural choices of statistics that are independent of the 
ordering of are 



T ^ 1 - - [(1 - A,y + (1 - A,y + (1 - A^f] . 

Both S and T statistics can be viewed as the suitably 
defined "distance" to the vertex {Ai,A2,A3) = (0,0,0). 
A third obvious choice, {Aj + Al + ^|)/3, is just 2S - T. 
To test alignment of the quadrupole and octopole planes 
with one another we quoted the S statistic numbers; T 
gives similar results. 

Alternatively, generalizing the definition in 22J, one 
can find, for each ^, the choice, n^, of z axis that maxi- 
mizes the angular momentum dispersion 



Li 



(13) 



One can then compare the maximized value with that 
from simulated isotropic skies [14]. Because £ = 2,3 are 
both planar (the quadrupole trivially so, the octopole be- 
cause the three planes of the octopole are nearly paral- 
lel) , the direction that maximizes the angular momentum 
dispersion of each is nearly the same as the (average) di- 
rection of that multipole's planes. Thus, the alignment 
of the octopole and quadrupole can be seen either from 
the S statistic, or by looking at the alignment of n2 with 

713. 

To test alignments of multipole planes with physical 
directions, we find the plane whose normal, n, has the 
largest dot product with the sum of the four quadrupole 
and octopole area vectors [T3]. Again, since • n is 
defined only up to a sign, we take the absolute value of 
each dot product. Therefore, we find the direction h that 



maximizes 



S 



^ Ni 



(14) 



1 



S=-{A, + A2 + As) 



and 



(12) 



C. Summary 

The study of alignments in the low-^ CMB has found 
a number of peculiarities. We have shown that the align- 
ment of the quadrupole and octopole planes is inconsis- 
tent with Gaussian, statistically isotropic skies at least 
at the 99% confidence level. Further a number of (possi- 
bly related) alignments occur at 95% confidence levels or 
greater. Put together these provide a strong indication 
that the full sky CMB WMAP maps are inconsistent with 
the standard cosmological model at large angles. Even 
more peculiar is the alignment of the quadrupole and oc- 
topole with solar system features (the ecliptic plane and 
the dipole). 

This is strongly suggestive of an unknown systematic in 
the data reduction; however, careful scrutiny has revealed 
no such systematic (except the mentioned modification 
of the radiometer gain model, that lead to a reduction 
of ecliptic alignment); see Sees. VC and VD for further 
discussion of the data analysis and instrumental explana- 
tions. We again stress that these results hold for full sky 
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maps; maps that are produced through combination of 
the individual frequency maps in such a way as to remove 
foregrounds. An alternative approach that removes the 
need for full sky maps is presented in the next section. 



IV. TWO-POINT ANGULAR CORRELATION 
FUNCTION 

The usual CMB analysis solely involves the spher- 
ical harmonic decomposition and the two-point angu- 
lar power spectrum. There are many reasons for this. 
Firstly, when working with a statistically isotropic uni- 
verse the angular power spectrum contains all of the 
physical information. Secondly, the standard theory pre- 
dicts the aim and their statistical properties, through the 
Ci, thus the spherical harmonic basis is a natural one to 
employ. Finally, as measured today the angular size of 
the horizon at the time of last scattering is approximately 
1 degree. Since 6'(deg) ~ 200/^ the causal physics at the 
surface of last scattering leaves its imprint on the CMB 
on small scales, 9 ^ 1° or £ > 100. The two-point angu- 
lar power spectrum focuses on these small scales, making 
it a good means of exploring the physics of the last scat- 
tering surface. The tremendous success of the standard 
model of cosmology has been the agreement of the theory 
and observations on these small scales, allowing for the 
precise determination of the fundamental cosmological 
parameters [3D]. 

The two-point angular correlation function provides 
another means of analyzing CMB observations and 
should not be ignored even if, in principle, it contains the 
same information as the angular power spectrum. Thus, 
even in the case of full sky observations and/or statisti- 
cal isotropy there are benefits in looking at the data in 
different ways. The situation is similar to a function in 
one dimension where it is widely appreciated that fea- 
tures easily found in the real space analysis can be very 
difficult to find in the Fourier transform, and vice versa. 
Furthermore, the two-point angular correlation function 
highlights behavior at large angles (small £); the opposite 
of the two-point angular power spectrum. Thus the angu- 
lar correlation function allows for easier study of the tem- 
perature fluctuation modes that are super-horizon sized 
at the time of last scattering. Finally, the angular corre- 
lation function in its simplest form is a direct pixel based 
measure (see below). Thus it does not rely on the recon- 
struction of contaminated regions of the sky to employ. 
This makes it a simple, robust measure even for partial 
sky coverage. 



A. Definition 

Care should be taken when discussing statistical quan- 
tities of the CMB and their estimators. It rarely is in the 
literature. Here we follow the notation of Copi et al. [5S], 



also see [3T] . The two point angular correlation function, 

C(ei,e2) = (T(ei)T(e2)), (15) 

is the ensemble average (represented by the angle brack- 
ets, (•)) of the product of the temperatures in the direc- 
tions ei and 62. Unfortunately we only have one universe 
to observe so this ensemble average cannot be calculated. 
Instead we average over the sky so that what we mean 
by the two-point angular correlation function is a sky av- 
erage, 

C{9) = T{i^)T{e2), (16) 

where the average is over all pairs of pixels with ei • 
62 = cos 9. This is a pixel based quantity and can be 
calculated for any region of the sky (of course not all 
separations 9 may be represented on a given patch of the 
sky, depending on its geometry). 



B. Missing angular power at large scales 

Spergel et al. [2j found that the two-point correlation 
function nearly vanishes on scales greater than about 60 
degrees, contrary to what the standard ACDM theory 
predicts, and in agreement with the same finding ob- 
tained from COBE data about a decade earlier [35]. 

We have revisited the angular two-point function in 
the 3-yr WMAP data in Ref. [55] and the 5-yr WMAP 
data in [3T]; see Fig. [s] From this figure we qualitatively 
see the following: 

• All of the cut-sky map curves are very similar to 
each other, and they are also very similar to the 
Legendre transform of the pseudo-Cf estimate of 
the angular power spectrum, which is not surpris- 
ing given that the two are formally equivalent [33) . 
Meanwhile the full-sky ILC C{9) and the Legen- 
dre transform of the maximum likelihood estimator 
(MLE) of the Cg agree well with each other, but not 
with any of the others. 

• The most striking feature of the cut-sky (and 
pseudo-Cf) C{9), is that all of them are very 
nearly zero above about 60°, except for some anti- 
correlation near 180°. This is also true for the full- 
sky curves, but less so. 

In order to be more quantitative about these obser- 
vations we adopt the 5*1/2 statistic introduced by the 
WMAP team [2] which quantifies the deviation of the 
two-point correlation function from zero, 

5i/2 = j^'\c{9)f d{cos9). (17) 

Spergel et al. [5] found that only 0.15% of the param- 
eter sets in their Markov chain of ACDM model CMB 
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FIG. 5. Two-point angular correlation function, C{9), com- 
puted in pixel space, for three different bands masked with 
the KQ75 mask (from WMAP 5 year data). Also shown is 
the correlation function for the ILC map with and without 
the mask, and the value expected for a statistically isotropic 
sky with best-fit ACDM cosmology together with 68% cosmic 
variance error bars. Even by eye, it is apparent that masked 
maps have C{9) that is consistent with zero at ^ > 60 deg. 
We also show the C{6) computed from the "official" pub- 
lished maximum likelihood estimator-based Ce- Clearly, the 
MLE-based Ct, as well as C{9) computed from the full-sky 
ILC maps, are in significant disagreement with the angular 
correlation function computed from cut-sky maps. Adopted 
from Ref. Ml- 



skies had lower values of Si /2 than the observed one- year 
WMAP sky. 

Applying this statistic we have found that the two- 
point function computed from the various cut-sky maps 
shows an even stronger lack of power, for WMAP 5 year 
maps significant at the 0.037%-0.025% level depending on 
the map used; see Fig. ([s]). However, we also found that, 
while C{9) computed in pixel space over the masked sky 
agrees with the harmonic space calculation that uses the 
pseudo-Cf estimator, it disagrees with the Ci obtained 
using the MLE (advocated in the 3rd year WMAP re- 
lease 15). The MLE-based Ce lead to C{9) that is low 
(according to the 5*1/2 statistic) only at the 4.6% level. 

There are actually two interesting questions one can 
ask: 

(i) Is the correlation function measured on the cut sky 
compatible with cut-sky expectation from the Gaus- 
sian random, isotropic underlying model? 

(ii) Is the reconstruction of the full sky correlation func- 
tion from partial information compatible with the 
expectation from the Gaussian random, isotropic 
underlying model? 

Our results refer to the first question above. The second 
question, while also extremely interesting, is more diffi- 
cult to be robustly resolved because the reconstruction 



FIG. 6. The two-point angular correlation function from the 
WMAP 5 year resuhs. Plotted are C{d) for the ILC calcu- 
lated separately on the part of the sky outside the KQ75 cut 
(dashed line), inside the KQ75 cut (dotted line), and on the 
part of the sky with at least on point inside the KQ75 cut 
(dotted-dashed line). For better comparison to the full-sky 
C{9) (solid line), the partial-sky C{9) have been scaled by the 
fraction of the sky over which they are calculated. Adopted 
from ReL [Pl- 



uses assumptions about statistical isotropy (see the next 
subsection) . 

The little large-angle correlation that does appear in 
the full sky maps (for example the solid, black line in 
Fig. [5| is associated with points inside the masked region. 
Shown in Fig. [6] are the normalized contributions to C{9) 
from different parts of the map. In particular, we see that 
almost all of the contribution to the full sky two-point 
angular correlation function comes from correlations with 
at least one point inside the masked region. Conversely, 
there is essentially no large-angle correlation for points 
outside the masked region and even very little among the 
points completely inside the mask. We also see that all 
the curves cross zero at nearly the same angle, 9 ~ 90°. 
We have no explanation for these results though they 
may point to systematics in the data. 



C. Alternative Statistics 

The two-point angular correlation function, C{9), as 
defined above in Eq. ( 15 ) is a simple pixel based measure 



of correlations. It makes no assumptions about the un- 
derlying theory, which can be taken as a feature or as a 
flaw. On the positive side, Eq. ( 15 ) does not assume that 



the standard model is correct and try to "force" it on the 
data. On the negative side we are not utilizing the full 
information available when comparing to the standard 
model. 

Various approaches have been taken to incorporate the 
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standard model in the analysis. For example, Hajian [33] 
defined a statistic that explicitly takes into account the 
covariance in the quantity C{9): 

/ c{e)F-^{e, e')C{e')d{cos e)d{co& 9'), 

(18) 

where F is the aforementioned covariance 

Fie,9') ^ {[Ci9) - {Ci9))] [C{9') - {Ci9'm, (19) 

and as usual the angle brackets denote an ensemble av- 
erage. Note that in the limit that C{9) is uncorrelated 
then A{l/2) = 5*1/2. Clearly this statistic relies on a 
model to calculate F{9,9'). With this statistic and as- 
suming the concordance model it is found that less than 
1% of realizations of the standard model have a A(0.53) 
less than those found for the masked skies. For the full 
sky ILC map approximately 8% of realizations have a 
smaller value. Though less constraining than making no 
assumptions about the theory through the use of the 5*1/2 
statistic, these results are consistent with those we pre- 
viously found. 

Another approach advocated by Efstathiou et al. [3S] 
is to reconstruct the full-sky C{9) from the partial sky 
and compare the reconstructed full-sky C{9) (using, say 
5*1/2) to the predictions of the model. This approach 
employs the usual map making algorithm on the low-^ 
spherical harmonic coefficients, aim- In this approach 
it is assumed that the statistical properties of the aim 
above some £max are known. In particular it assumes 
there are no correlations between the aim with i < ^,nax 
and those with i > ^max- The method is very similar to 
a maximum likelihood analysis (see I36j . for example). 
As we have seen above (e.g. in Fig. pi it is not surpris- 
ing that this approach will be consistent with the full-sky 
ILC results, as these authors have verified. As mentioned 
in the previous subsection, however, that this procedure 
poses a different question of the data than has been ad- 
dressed by the 5*i/2 statistic applied to a masked sky. 
With the map-making technique a full-sky map is con- 
structed that is consistent with the sky outside the mask 
but relies on assumptions to fill in the masked region. 
As is clear from Fig. [6] it is precisely the region inside 
the mask that is introducing correlations with the region 
outside the mask. Thus, the assumptions required to al- 
low filling in the masked region also produce two-point 
angular correlations. 

Whether or not reconstructing the full sky is a "more 
optimal" approach than direct calculation of the cut-sky 
C{9) is moot, but likely depends on the actual (rather 
than the assumed) statistical properties of the underly- 
ing fluctuations, as well as on the particular realization 
of those distributions. What is important is to make 
an "apples to apples" comparison between the observed 
sky and simulated realizations of the ensemble of possi- 
ble skies. As we have shown, the pixel-based two-point 
correlation function on the region of the sky outside a 
conservative galactic mask is inconsistent with the pre- 
dictions of the standard ACDM model for the identical 



pixel-based two-point correlation function on the iden- 
tically masked sky. The fact that the full-sky analysis 
shows less statistical significance is not in contradiction 
with the cut-sky result, although it may eventually help 
in pointing to a cause of this anomaly. 



D. Summary 

The striking feature of the two-point angular correla- 
tion function as seen in Fig. [5] is not that it disagrees 
with ACDM (though it does at > 90% C.L.) but that at 
large angles it is nearly zero. This lack of large angle cor- 
relations is unexpected in inflationary models. The 5i/2 
statistic quantifies the deviation from zero and shows a 
discrepancy exists at more than 99.9% C.L. Equally strik- 
ing is the fact that the little correlation that does exist 
in the full sky ILC map, or equally in the MLE esti- 
mated aim, comes from correlations between the masked 
foreground region and the expectedly cleaner CMB re- 
gions of the sky. Thus this residual correlation, which 
still is discrepant with generic inflationary predictions at 
about 95% C.L., comes from the reconstruction proce- 
dure. This surprising lack of large angle correlation out- 
side the masked region remains an open problem. 

We also note that the vanishing of power is much more 
apparent in real space (as in C{9)) than in multipole space 
(as in Ci). The harmonic-space quadrupole and octopole 
are only moderately low (e.g. [37]), and it is really a range 
of low multipoles that conspire to make up the vanishing 
C{9). Specifically, as discussed in [31 , there is a cancella- 
tion between the combined contributions of C2,---,G^ and 
the contributions of Ci with £ > 6. It is this conspiracy 
that is most disturbing, since it violates the independence 
of the Ci of different I that defines statistical isotropy. 

In [3T] we therefore explored the possibility that the 
vanishing of C{9) could be explained simply by changing 
the values of the theoretical \ow-l Ci, as might be the 
result, say, of a modified inflaton potential. In particular, 
we replaced C2 through C20 in the best-fit ACDM model 
with the values extracted from the cut-sky ILC five-year 
map. From these Ci, 200, 000 random maps were created, 
masked, and 5i/2 computed. Under the assumptions of 
Gaussianity and statistical isotropy of these Ci only 3 
per cent of the generated maps had 5i/2 less than the 
cut-sky ILC5 value. Thus, even if the Ci are set to the 
specific values that produce such a low 5*i/2, a Gaussian 
random, statistically isotropic realization is unlikely to 
produce the observed lack of large angle correlations at 
the 97% C.L. Moreover, in work in progress we show that 
almost all of those 3% are skies with several anomalously 
low Ci — not at all the sky we see. Only a tiny fraction 
of the 3% represent skies in which most of the individual 
Ci were close to the ACDM prediction but they conspire 
to cancel one another in the large angle C{9). This shows 
that either (i) the low-£ Ce are correlated, contrary to the 
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assumption of statistical isotrop}|jor (ii) our Universe is 
an extremely unlikely realization of whatever statistically 
isotropic model one devises. 

It is for this reason that theoretical efforts to explain 
"low power on large scales" must focus on explaining the 
low C{9) at ^ 60 deg, rather than the low quadrupole. 

Finally, one might ask if the observed lack of correla- 
tion and the alignment of quadrupole and octopole are 
correlated. This issue was studied by Rakic & Schwarz 
[55] for the full sky and by Sarkar et al. [31] for the cut 
sky case. In both cases, it was shown that low power 
and alignments are uncorrelated — i.e. that having one 
does not imply a larger or smaller probability of hav- 
ing the other. This was shown by applying a Monte- 
Carlo analysis to sky realizations with the underlying 
standard Gaussian random, statistically isotropic cosmo- 
logical model, without any further constraints. Thus one 
might view the 99.6% C.L. of quadrupole-octopole align- 
ment presented in the previous section and the 95% C.L. 
for lack of correlation in full-sky maps reported in this 
section as statistically independent. 



V. QUEST FOR AN EXPLANATION 

Understanding the origin of CMB anomalies is clearly 
important. Both the observed alignments of the \ow-£ 
full-sky multipoles, and the absence of large-angle corre- 
lations (especially on the galaxy-cut sky) are severely in- 
consistent with predictions of standard cosmological the- 
ory. There are four classes of possible explanations: 

1. Astrophysical foregrounds, 

2. Artifacts of faulty data analysis, 

3. Instrumental systematics, 

4. Theoretical/cosmological. 

In this section, we review these four classes of expla- 
nations, giving examples from each. First, however, we 
discuss two generic ways to break statistical isotropy and 
affect the intrinsic (true) CMB signal — additive and 
multiplicative modulations — and illustrate in general 
terms why it has been so difficult to explain the anoma- 
lies. 



A. Additive vs. multiplicative effects 

Why is it difficult to explaining the observed CMB 
anomalies? There are three basic reasons: 



Though this appears to be an unhkely explanation since corre- 
lations between the Ci generically increase the variance of the 
Si/2 statistic |33) . 



• Most explanations work by adding power to the 
large-angle CMB, while the observed anisotropics 
actually have less large-scale power, and particu- 
larly less large angle correlation, than the ACDM 
cosmological model predicts. 

• Unaccounted for sources of CMB fluctuations in 
the foreground, even if possessing/causing aligned 
low-^ multipoles of their own, cannot bring un- 
aligned statistically isotropic cosmological pertur- 
bations into alignment. Therefore, aligned fore- 
grounds as an explanation for alignment work only 
if the cosmological signal is subdominant, thus ex- 
acerbating the lack of large-angle correlations. 

• The alignments of the quadrupole and octopole are 
with respect to the ecliptic plane and near the 
dipole direction. It is generally difficult to have 
these directions naturally be picked out by any class 
of explanations (though there are exceptions to this 
— see the instrumental example below). 

Gordon et al. [40 , Rakic & Schwarz [38] and Bunn & 
Bourdon [31] explored generic "additive models" where 
the temperature modification that causes the alignment 
is added to the intrinsic temperature 

-^observed (^) -^intrinsic (e) + radd(e). (20) 

Here Tadd(e) is the additive term — perhaps contamina- 
tion by a foreground, perhaps an additive instrumental 
or cosmological effect. They showed that additive mod- 
ulations of the CMB sky that ameliorate the alignment 
problems tend to worsen the overall likelihood at large 
scales (though they may pick up offsetting positive likeli- 
hood contribution from higher multipoles) . The intuitive 
reason for this is that there are two penalties incurred by 
the additive modulation. First, since the power spectrum 
at low i is lower than expected, one typically needs to ar- 
range for an accidental cancellation between Ti„trinsic and 
Tadd; the cancellation moreover must leave an aligned 
quadrupole and octopole even though the quadrupole 
and octopole of Tlntrinsic are not aligned. (Very simi- 
lar reasoning argues against additive explanations of the 
suppression of large angle correlations.) Second, the sim- 
plest models for the additive contribution that are based 
on an azimuthally symmetric modulation of a gradient 
field can only affect m = multipoles around the pre- 
ferred axis, while as we mentioned earlier the observed 
quadrupole and octopole as seen in the preferred (dipole) 
frame are dominated by the m = £ components. 

In contrast to the additive models, the multiplicative 
mechanisms, where the intrinsic temperature is multi- 
plied by a spatially varying modulation, are phenomeno- 
logically more promising. As a proof of principle, a toy- 
model modulation [iO] 

7'obsorvod(e) =./[!+ W2>20(e)] T'intrinsic(e), (21) 

(where the modulation is a pure I20 along the dipole 
axis), can increase the likelihood of the WMAP data by 
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FIG. 7. A realization of the multiplicative model where the 
quadrupole (left column) and octopole (right column) exhibit 
an alignment similar to WMAP. First row: intrinsic (unmodu- 
lated) sky from a Gaussian random isotropic realization. Sec- 
ond row (single column): the quadrupolar modulation with 
/ = —1 and W2 — —7 (see Eq. (21l) in the dipole direction. 
Third row: the modulated sky of the observed CMB. Fourth 
row: WMAP full-sky quadrupole and octopole. Adopted from 
Ref. [ID]. 



a factor of exp (16/2) and, at the same time, increase the 
probability of obtaining a sky with more alignment (e.g. 
higher angular momentum statistic) 200 times, to 45%; 
see Fig. [rj Indeed, Groeneboom et al. [32], building on 
the work of Groeneboom & Eriksen [13] and Hanson & 
Lewis [44] and motivated by a model due to Ackerman 
et al. [15] , have identified a 9cr quadrupolar power asym- 
metry, recently confirmed by the WMAP team [IS]; this 
anomaly can, however, be fully explained by accounting 
for asymmetric beams [47] . Recently, Hoftuft et al. ^48j 
found a greater than 3-cr evidence for nonzero dipolar 
modulation of the power. 



B. Astrophysical explanations 

One fairly obvious possibility is that there is a perni- 
cious foreground that contaminates the primordial CMB 
and leads to the observed anomalies. Such foregrounds 
are, of course, additive mechanisms, in the sense of the 
preceding section, and so suffer from the shortcomings 
described therein. Moreover, most such foregrounds are 
Galactic, while the observed alignments are with respect 
to the ecliptic plane. One would expect that Galactic 
foregrounds should lead to Galactic and not ecliptic fore- 
grounds. This simple expectation was confirmed in |14| . 
where we showed that, by artificially adding a large ad- 
mixture of Galactic foregrounds to WMAP CMB maps, 
the quadrupole vectors move near the z-axis and the nor- 
mal into the Galactic plane, while for the octopole all 



three normals become close to the Galactic disk at 90° 
from the Galactic center. Therefore, as expected Galac- 
tic foregrounds lead to Galactic, and not ecliptic, corre- 
lations of the quadrupole and octopole (see also studies 
by [33 [SD]). 

Moreover, in |14| . we have shown that the known 
Galactic foregrounds possess a multipole vector struc- 
ture very different from that of the observed quadrupole 
and octopole. The quadrupole is nearly pure Y22 in the 
frame where the z-axis is parallel to the dipole (or li^^^.s) 
or any nearly equivalent direction), while the octopole is 
dominantly 1^33 in the same frame. Mechanisms which 
produce an alteration of the microwave signal from a rel- 
atively small patch of sky — and all of the recent pro- 
posals fall into this class — are most likely to produce 
aligned I20 and Y30. This is essentially because the low-f 
multipole vectors will all be parallel to each other, lead- 
ing to a Yfo in this frame. 

A number of authors have attempted to explain the 
observed quadrupole-octopole correlations in terms of a 
new foreground — for example the Rees-Sciama effect 
[FT] , interstellar dust [52], local voids or the 

Sunyaev-Zeldovich effect [54]. Most if not all of these 
proposals have a difficult time explaining the anomalies 
without severe fine tuning. For example, Vale |55j clev- 
erly suggested that the moving lens effect, with the Great 
Attractor as a source, might be responsible for the extra 
anisotropy; however Cooray & Seto fSS] have argued that 
the lensing effect is far too small and requires too large 
a mass of the Attractor. 

It is also interesting to ask if any known or unknown 
Solar system physics could lead to the observed align- 
ments. Dikarev et al. [53 |5H] studied the question of 
whether solar system dust could give rise to sizable levels 
of microwave emission or absorption. Surprisingly, very 
little is known about dust grains of mm to cm size in 
the Solar system, and their absorption/emission prop- 
erties strongly depend on their chemical composition. 
While iron and ice particles can definitely be excluded 
to contribute at significant levels, carboneous and sili- 
cate dust grains might contribute up to a few /iK close 
to the ecliptic plane, e.g. due to the trans-Neptunian ob- 
ject belt. Such an extra contribution along the eclip- 
tic could give rise to CMB structures aligned with the 
ecliptic, but those would look very different from the ob- 
served ones. On top of that, Solar system dust would 
be a new additive foreground and could not explain the 
lack of large angle correlations. Thus it seems unlikely 
that Solar system dust grains cause the reported large an- 
gle anomalies, nevertheless they are sources of microwave 
absorption and emission and may become important to 
precision measurements in the future. 

Finally, it has often been suggested to some of us in 
private communications that the anomalies may not re- 
fiect an unknown foreground that has been neglected, 
but rather the "mis-subtraction" of a known foreground. 
However, it has never quite been clear to us how this 
leads to the observed alignments or lack or large angle 
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correlations, and we are unaware of any literature that 
realizes this suggestion successfully. 



C. Data analysis explanations 

Most of the results discussed so far have been obtained 
using reconstructed full-sky maps of the WMAP obser- 
vations [H 122 EH]. In the presence of the sky cut of 
even just a few degrees, the errors in the reconstructed 
anisotropy pattern, and the directions of multipole vec- 
tors, are too large to allow drawing quantitative conclu- 
sions about the observed alignments [15]. These large 
errors are expected: while the power in the CMB (rep- 
resented, say, by the angular power spectrum Cg) can be 
accurately recovered since there are 2£+l modes available 
for each £, there are only 2 modes available for each mul- 
tipole vector; hence the cut-sky reconstruction is noisier. 
However the cut-sky alignment probabilities, while very 
uncertain, are consistent with the full-sky values [TH ISO] ; 
more generally, the alignments appear to be rather robust 
to Galactic cuts and foreground contamination [60j . 

A different kind of explanation of missing large-scale 
power, or missing large angle correlations, has been taken 
by Efstathiou et al. [35j who argued that maximum like- 
lihood estimators can be applied to the cut-sky maps to 
reliably and optimally reconstruct the CMB anisotropy 
of the whole sky; for a recent work that extends these 
ideas, see |33j . This approach yields more two-point cor- 
relation on large scales (S'1/2 ^ 8000 (/xif)'*, which is 
~ 5% likely) than the direct cut-sky pixel-based calcula- 
tion which gives S'1/2 ~ 1000 (l^K)'^ result and is ^ 0.03% 
likely. These authors then argue that the extremely low 
Si/2 obtained by the pixel-based approach applied to the 
cut sky is essentially a fluke, and the more reliable result 
comes from their maximum-likelihood reconstruction of 
the full sky. It may indeed be true that the pixel-based 
calculation is a suboptimal estimate of the full-sky C{9) 
for a statistically isotropic cosmology. However, quanti- 
ties calculated on the cut sky are clearly insensitive to 
assumptions about what lies behind the cut. We can 
only observe reliably the ~ 75% of the sky that was 
not masked, and that is where the large-angle two-point- 
correlation is near- vanishing. Any attempt to reconstruct 
the full sky must make assumptions about the statistical 
properties of the CMB sky, and would clearly be affected 
by the coupling of small-scale and large-scale modes — 
exactly what is necessary to have a sky in which S'1/2 is 
anomalously low, while the C^s are individually approx- 
imately consistent with the standard cosmology. 



away from the ecliptic more than those near the eclip- 
tic. While the corresponding variations in the noise per 
pixel are well known (e.g. as the number of observations 
per pixel, A'obs; see Fig. 3 in [T]), and its effects on the 
large-scale anomalies are ostensibly small, they could, in 
principle, be amplified and create the observed ecliptic 
anomalies. However a successful proposal for such an 
amplification has not yet been put forward. 

Another possibility is that an imperfect instrument 
couples with dominant signals from the sky to create 
anomalies. Let us review an example given in Ref. |40) : 
suppose that the instrumental response Tinstr(e) to the 
true sky signal T{e) is nonlinear 
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(22) 



Here / is an arbitrary normalization scale for the non- 
linearity of the response, and ai are arbitrary coefficients 
with ai = 1. If aiyi =/= then Ti^sti 7^ T and the 
observed temperature is a non-linear modulation of the 
true temperature. The dominant temperature signal for 
a differencing experiment such as WMAP is the dipole 
arising from our peculiar motion, T{e) = Tdip cos 9, with 
Tdip = 3.35mK and 6 the polar angle in the dipole frame. 
Taking / = Tdip, 
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where Pe are the Legendre polynomials. Note that with 
Q2 ^ Q!3 ^ 0(10^^), the 10^"^ dipole anisotropy is mod- 
ulated into a 10^^ quadrupole and octopole anisotropy 
which are aligned in the dipole frame with the m = 
multipole structure. Unfortunately (or fortunately!), 
WMAP detectors are known to be linear to much better 
than 1%, so this particular realization of the instrumen- 
tal explanation does not work. As an aside, note that 
this type of explanation needs to assume that the higher 
multipoles are not aligned with the dipole/ecliptic, and 
moreover, requires essentially no intrinsic power at large 
scales (that is, even less than what is observed). 

To summarize: even though the ecliptic alignments 
(and the north-south power asymmetry) hint at a sys- 
tematic effect due to some kind of coupling of an obser- 
vational strategy and the instrument, to date no plausible 
proposal of this sort has been put forth. 



D. Instrumental explanations 

Are instrumental artifacts the cause of the observed 
alignments (and/or the low large-scale power)? One pos- 
sible scenario would go as follows. WMAP avoids mak- 
ing observations near the Sun, therefore covering regions 



E. Cosmological Explanations 

The most exciting possibility is that the observed 
anomalies have primordial origin, and potentially inform 
us about the conditions in the early universe. One ex- 
pects that in this case the alignments with the dipole, or 
with the solar system, would be statistical fiukes. 
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The breaking of statistical isotropy implies that the 
usual relation {a}^aim) = CiSwSmm' does not hold any 
more; instead 



dark matter perturbations 



= a 



ii'nim' ! 



(24) 



where the detailed form of the quantity on the right-hand 
side is model-dependent (see the more detailed discus- 
sions in [STHB3]). 

There are many possibilities for how the absence of 
statistical isotropy might arise. For example, in a non- 
trivial spatial topology, the fundamental domain would 
not be rotationally invariant, and so the spherical har- 
monics (times an appropriate radial function) would not 
be a basis of independent eigenmodes of the fluctuations. 
This would certainly lead to a correlation among aim of 
different m and also £, although not necessarily to aligned 
multipoles. The hope would be that the shape of the fun- 
damental domain would lead to these alignments, while 
a lower density of states at long wavelength (compared 
to the covering space) would lead to the absence of large 
angle correlations. No specific model has been suggested 
to accomplish all of those, and the matter is complicated 
by known bounds on cosmic topology which force 

the fundamental domain to be relatively large. 

Alternately, in the early universe, asymmetry in the 
stress-energy tensor of dark energy [5^], or a long- 
wavelength dark energy density mode with a gradient 
in the desired direction |40j, could both imprint the ob- 
served alignments via the integrated Sachs- Wolfe mech- 
anism; but it is hard to see how they would explain the 
lack of large angle correlations. The authors of [70 have 
put forward a model where the Sachs- Wolfe contribution 
to \ow-£ multipoles is partly cancelled by the Integrated 
Sachs- Wolfe contribution, but which still fails to explain 
the lowness of S'1/2 or the alignments of low-^ multipoles. 
While models where the anomalies are caused by break- 
ing of the statistical isotropy [TTl [72] are especially well 
studied, see [77 for the equally interesting possibility 
that the translational invariance is broken. 

A commonly used mechanism to explain such anoma- 
lies are inflationary models that contain implicit break- 
ing of isotropy [IHl ITiUSOj . Pontzen [81 helpfully shows 
temperature and polarization patterns caused by various 
classes of Bianchi models that explicitly break the statis- 
tical isotropy. However, outside of explaining the anoma- 
lies, the motivation for these anisotropic models is not 
compelling and they seem somewhat contrived. More- 
over, the authors have not investigated whether the low 
Si/2 is also observed. Nevertheless, given the large-scale 
CMB observations, as well as the lack of fundamental 
theory that would explain inflation, investigating such 
models is well worthwhile. 

A very reasonable approach is to describe breaking of 
the isotropy with a phenomenological model, measure the 
parameters of the model, and then try to draw inferences 
about the underlying physical mechanism. For exam- 
ple, a convenient approach is to describe the breaking of 
isotropy via the direction-dependent power spectrum of 



P{k) = A(k) 



im 



(25) 



where k — \k\, gim{k) quantifies the departure from sta- 
tistical isotropy as a function of wavenumber fc, and A(k) 
refers to the statistically isotropic part of the power spec- 
trum. In this model, the power spectrum, normally con- 
sidered to depend only on scale fc, now depends on direc- 
tion in a parametric way. Statistically significant finding 
that gim 7^ for any (£, m) would signal a violation of 
statistical isotropy. 

As with the other attempts to explain the anomalies, 
we conclude that, while there have been some interesting 
and even promising suggestions, no cosmological expla- 
nation to date has been compelling. 



F. Alignment explanations: what next 

While future WMAP data is not expected to change 
any of the observed results, our understanding and analy- 
sis techniques are likely to improve. The most interesting 
test will come from the Planck satellite, whose temper- 
ature maps, obtained with a completely different instru- 
ment and observational technique than WMAP, could 
shed significant new light on the alignments. Moreover, 
polarization information could be extremely useful in dis- 
tinguishing between different models and classes of ex- 
planations in general; for example, Dvorkin et al. |83) ex- 
plicitly show how polarization information expected from 
Planck can help identify the cause of the alignments. Fi- 
nally, one could use the large-scale structure (i.e. galaxy 
distribution) data on the largest observable scales from 
surveys such as Dark Energy Survey (DES) and Large 
Synoptic Survey Telescope (LSST) to test cosmological 
explanations; see e.g. [Ml 185) . 



VI. EXPLANATIONS FROM THE WMAP 
TEAM 

In their seven year data release the WMAP team ex- 
plicitly discusses several CMB anomalies [15] including 
the two main ones described in this review. For the first 
major issue — the alignment of low multipoles with each 
other — the WMAP team agrees that the alignment is 
observed and argue, based on work by Francis and Pea- 
cock [86j, that the integrated Sachs- Wolfe (ISW) contri- 
bution of structures at small redshifts (z <C 1) could be 
held responsible. There are serious problems with this 
argument. Firstly, the ordinary Sachs- Wolfe (SW) effect 
typically dominates at these £ over the ISW. Thus, only 
if the ordinary SW effect on the last scattering surface 
is anomalously low will the ISW contribution dominate. 
Secondly, though the ISW may lead to alignment of the 
quadrupole and octopole it is not an explanation for the 
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observed Solar system alignments. This alignment would 
need to be an additional statistical fluke. Finally, this ex- 
planation does nothing whatsoever to mitigate the lack 
of large scale angular correlation because the ISW effect 
acts as an additive component and should be statisti- 
cally uncorrelated from the primordial CMB. Therefore, 
even if the ISW reconstruction is taken as reliable, this 
argument would imply; 

1. an accidental downward fluctuation of the SW suf- 
flcient for the ISW of local structure to dominate 
and cause an alignment, and 

2. an accidental cancellation in angular correlation be- 
tween the SW and ISW temperature patterns. 

Neither the WMAP team nor Francis and Peacock esti- 
mate the likelihood of these two newly created puzzles. 

Regarding the second major issue — the lack of an- 
gular correlation — the WMAP team refers to a recent 
work by Efstathiou, Ma and Hanson [5S] who argue that 
quadratic estimators are better estimates of the full sky 
from cut-sky data and are in better agreement with the 
concordance model. While these estimators have been 
shown to be optimal under the assumption of statistical 
isotropy, it is unclear why they should be employed when 
this assumption is to be tested. (For a contrary view, 
see [33].) The pixel based estimator applied to the cut 
sky in our analysis does not rely on statistical isotropy 
and it is more conservative as it does not try to recon- 
struct temperature anisotropics inside the cut. Finally, 
our claims that the pixel based cut-sky two-point corre- 
lation function is highly anomalous rely on comparisons 
to the identical correlation function calculated on sim- 
ulated cut skies. Whether or not it is a good estimate 
of the full-sky correlation function is answering a differ- 
ent question, (see Sec. IV B). Conversely, the estimator 
suggested in [35] assumes that whatever is within the cut 
can be reconstructed reliably by truncating the number of 
multipole moments considered. The latter logic is equiv- 
alent to the assumption of the statistical independence of 
low and high multipoles, which is exactly a consequence 
of statistical isotropy. 

These arguments from the WMAP team offer neither 
new nor convincing explanations of the observed anoma- 
lies discussed in this review. At best they replace one set 
of anomalies for another. 



VII. CONCLUSIONS 

The CMB is widely regarded as offering strong sub- 
stantiating evidence for the concordance model of cos- 
mology. Indeed the agreement between theory and data 
is remarkable — the patterns in the two point correla- 
tion functions (TT, TE and EE) of Doppler peaks and 
troughs are reproduced in detail by fitting with only six 
(or so) cosmological parameters. This agreement should 
not be taken lightly; it shows our precise understanding 



of the causal physics on the last scattering surface. Even 
so, the cosmological model we arrive at is baroque, re- 
quiring the introduction at different scales and epochs 
of three sources of energy density that are only detected 
gravitationally — dark matter, dark energy and the in- 
flaton. This alone should encourage us to continuously 
challenge the model and probe the observations particu- 
larly on scales larger than the horizon at the time of last 
scattering. 

At the very least, probes of the large-angle {\ow-£) 
properties of the CMB reveal that we do not live in a 
typical realization of the concordance model of inflation- 
ary ACDM. We have reviewed a number of the ways in 
which that is true: the peculiar geometry of the £ = 2 
and 3 multipoles — their planarity, their mutual align- 
ment, their alignment perpendicular to the ecliptic and 
to the dipole; the north-south asymmetry; and the near 
absence of two-point correlations for points separated by 
more than 60°. 

If indeed the observed £ = 2 and 3 CMB fluctuations 
are not cosmological, one must reconsider all CMB re- 
sults that rely on the low £, e.g. the measurement of 
the optical depth from CMB polarization at low £ or the 
spectral index of scalar perturbations and its running. 
Moreover, the CMB-galaxy cross-correlation, which has 
been used to provide evidence for the Integrated Sachs- 
Wolfe effect and hence the existence of dark energy, also 
gets contributions from the lowest multipoles (though 
the main contribution comes from slightly smaller scales, 
i ^ 10). Indeed, it is quite possible that the underlying 
physical mechanism does not cut off abruptly at the oc- 
topole, but rather affects the higher multipoles. Indeed, 
several pieces of evidence have been presented for anoma- 
lies at Z > 3 (e.g. [53 HE])- One of these is the parity of 
the microwave sky. While the observational fact that 
the octopole is larger than the quadrupole (C3 > C2) is 
not remarkable on its own, including higher multipoles 
(up to ^ ~ 20) the microwave sky appears to be par- 
ity odd at a statistically significant level (since WMAP 
5yr) [891 - I9T] . It is hard to imagine a cosmological ex- 
planation for a parity odd universe, but the same holds 
true for unidentified systematics or unaccounted astro- 
physical foregrounds, especially as this recently noticed 
puzzle shows up in the very well studied angular power 
spectrum. 

While the further WMAP data is not expected to 
change any of the observed results, our understanding 
and analysis techniques are likely to improve. Much 
work remains to study the large-scale correlations using 
improved foreground treatment, accounting even for the 
subtle systematics, and in particular studying the time- 
ordered data from the spacecraft. The Planck experi- 
ment will be of great importance, as it will provide maps 
of the largest scales obtained using a very different ex- 
perimental approach than WMAP — measuring the ab- 
solute temperature rather than temperature differences. 
Polarization maps, when available at high enough signal- 
to- noise at large scales (which may not be soon) , will be a 
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fantastic independent test of the alignments, as each ex- 
planation for the alignments, in principle, also predicts 
the statistics of the polarization pattern on the sky. 
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